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ABSTRACT 

We study the effect of active galactic nuclei (AGN) heating on the chemodynami- 
cal evolution of elliptical galaxies and their X-ray and optical properties using high- 
resolution A-dominatcd cold dark matter cosmological simulations. Our model con- 
siders an AGN as being "active" when a convergent gas inflow condition exists within 
a galaxy's nucleus; otherwise, the AGN is assumed to remain dormant. This induces 
a self-regulated activity for the AGN, the result of which leads to a stable hot corona 
and the suppression of significant late-time star formation - characteristics not en- 
countered in traditional chemodynamical models of ellipticals. These properties of our 
AGN heating model leads to a system consistent with both the X-ray and optical 
properties of comparable elliptical galaxies. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: formation — galaxies: 
evolution — galaxies: stellar content 



1 INTRODUCTION 

, One of the primary science drivers for the next generation 
' of ground- and space-based observational facilities which 
spread over multi-wavelength is an understanding of the 
, physics of galaxy formation and evolution. In the case of el- 
' liptical galaxies, the optical regime provides constraints on 
\ the properties of the underlying stellar populations, while 
. the X-ray regime yields insights into the physical conditions 
' of t he associated hot in t erstella r (coronal) medium. 

iKawata fc Gibson I J2003bl KG03b) presented a first 
attempt at explaining both the X-ray and optical prop- 
erties of ellipticals via the use of self-consistent cosmo- 
logical simulations. Using a standard "recipe" for galaxy 
formation, they found that radiative cooling is important 
to interpret the observed X-ray luminosity (Lx), temper- 
ature (Tx), and metallicity ([Fe/H]x) of the hot gas of 
elliptical galaxies. However, these models were subject to 
an unavoidable serious problem, in that the cooled gas 
unavoidably le d to excessive star formation at low red- 
shift (see also 'Suginohara & Ostriker ' 'l99^; iLewi^e^Lj 
[2000 ; Kav Thom as fc Thcuns 2003; T ornator e et al. l2003t 
iMeza et al. ll20o3:lBorgani et al. II2CIO4I1 . This late-time star 
formation led to associated stellar populations which were 
too blue with respect to observations. 

KG03b suggested that a heating mechanism not in- 
cluded in their simulations would be required to suppress 
this enhanced cooling and consequent star formation. In 
this paper, we examine one such potential mechanism - 
heating by an active galactic nucleus (AGN). In such a 
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picture, one might expect an AGN could become "active" 
only when being fed "fuel" - i.e., when cold gas inflows 
into the central region, the AGN would activate and heat 
the surrounding gas, and potentiall y balance the associ- 
ated radiative cooling (e.g. [ Churazov ^^^^jboO^ jNulsen I 
12004 iGranato et al. 1120041: IScannapieco fc Ohll2004^ . Once 
the infalling fuel source is quenched by the AGN heating, 
the AGN reverts to its quiescent state, and cycle can be- 
gin anew. This self-regulation provides heating (suppress- 
ing star formation) without compromising cooling (required 
to recover the empirical X-ray properties). Recent high- 
resolution X-ray images taken by Chandra reveal that the 
central hot gas of ellipticals is not smoothly distributed, but 
possesses cavities on scales comparable to the radio cmis- 
sion (e.g.feohringer et al. ,.19 93: Finogucnov & Jones 20o3; 
I Jones et al. 120021) . The X-rav holes are also seen on the clus- 
ter scale (> 50 k pc), and are often coincident with extended 
radio lobes (e.g. |McNMnBxa^.et_aL I I2OOOI : lAllen et al "I I2OOII : 
iBlanton et al. l200ll : lFuiita et al. l2002l : lHeintz et al. I200I1 . 
These features are considered to be a relic of the influence 
of AGN activity upon the hot gas in ellipticals and clusters 
JChurazov et al. 11200 ih . 

Our current study examines the effect of AGN heating 
on both the optical and X-ray properties of giant elliptical 
galaxies in the context of a self-consistent treatment 
of cosmological evolution. Several earlier studies of the 
effect of AGN heating on the h ot ga s of elliptic als can 
be se e n in [B inney ^^TaborJ il995(); 'Churazov et al. 
i2001 ); ICiotti fc Ostriker I J2001I): Quilis, Bower fc Balogl^ 
12001 ): iBrighenti fc Mathews I (l2002i'): .Bruggen fc KaiseT 
120021): ISaxton Bicknell fc Sutherland! ~ |2002[): 

iBasson fc Alexander I ll2003l) : lOmma et all J2004) : 
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Figure 1. Comparison of tlie simulated and observed (crosses 
with error bars) Lx — Tx relations (upper) and [Fe/H]x — Lx 
relations (lower). The square/circle/triangle indicates the predic- 
tions of Model 1/2/2L. Open stars in the lower panel show the 
predicted mean metallicity of the stellar component for the two 
models. The observational data are from MOMOO. The shaded 
area represents the range of Lx and Tx for a sample of elliptical 
galaxies with low-Lx and high-Tx (group 5 in MOMOO), and 90 
% confidence limit of [Fe/H]x from the composite X-ray spec- 
trum of these galaxies. The solid star corresponds to one of these 
galaxies, NGC 3923. 



iHoeft fc Briieggenl l)2004) . While all important con- 
tributions to the field, these previous studies were (a) 
restricted in their temporal coverage to <15% of a Hubble 
time, and (b) did not trace star formation or mass accretion 
self-consistently within a A-dominated cold dark matter 
(ACDM) cosmology. Our cosmological simulation makes 
it possible to investigate whether or not AGN heating 
can lead to (a) a stable condition for the hot gas, and (b) 
suppress late-time star formation. 



2 METHODS 

Our simulations wer e carried out using the gal actic chemo- 
dynamics code GCD+ jKawata fc Gibson Il2003 a|) with initial 
conditions identical to those described in KG03b. GCD+ is a 
three-dimensional tree A''-body /smoothed particle hydrody- 
namics (SPH) code which incorporates self-gravity, hydrody- 
namics, radiative cooling, star formation, supernovae (SNe) 
feedback, and metal enrichment. GCD+ takes account of the 
chemical enrichment by both Type II (SNe II) and Type la 
(SNe la) SNe, mass-loss from intermediate mass stars, and 
follows the chemical enrichment history of both the stellar 



and gas components of the system. Here, we briefly describe 
our cosmological simulation models. 

We adopt a ACDM cosmology (r2o=0.3, Ao=0.7, 
S7b=0.019/i^^, h—0.7, and as=0.9) and use a multi- 
resolution technique to achieve high-resolution in the re- 
gions of interest, including the tidal forces from neighbour- 
ing large-scale structures. The initial conditions for the sim- 
ulations arc constructed using the public software GRAFIC2 
^Bertsc hinger 2001). Gas dynamics and star formation are 
included only within the relevant high-resolution region 
(~12 Mpc at z—Q); the surrounding low- resolution region 
(~43 Mpc) contributes to the high-resolution region only 
through gravity. Consequently, the initial condition consists 
of total 190093 dark matter particles and 134336 gas parti- 
cles. The mass and softening length of individual gas (dark 
matter) particles in the high-resolution region are 5.86 x 10^ 
(3.95 X 10*) Mq and 2.27 (4.29) kpc, respectively. KG03b 
found an appropriate elliptical galaxy analog in the high- 
resolution region, which acts as the focus for this study. The 
total virial mass of this galaxy is 2 x 10^^ M©. The target 
galaxy is relatively isolated, with only a few low-mass satel- 
lites remaining at z=0. Fig. 1 of KG03b shows the morpho- 
logical evolution of dark matter in the simulation volume 
and the evolution of the stellar component of the target 
galaxy. The galaxy forms through conventional hierarchical 
clustering between redshifts 2=3 and 2=1. The morphology 
has not changed dramatically since 2=1. 

KG03b presented the results of three different radiative 
cooling and SNe feedback models. We use one of their models 
as a reference model for comparison with this study; Model 
1 here corresponds to Model C of KG03b, which includes 
cooling and strong SNe feedback.^ An additional model 
(Model 2) is co nstructed, to study the effect of heating by an 
AGN (see also lSoringel. Pi matteo. fc Hernauist 112004311^) . 
In Model 2, strong thermal SNe feedback (as in Model 1) 
is also adopted. In addition, the most bounded star particle 
in the target galax;y at 2 = 1 is assumed to be the heating 
source - i.e., the AGN "particle".'^ Since any AGN heating 
necessarily requires a fuel (gas) source, we assume that the 
AGN heating is active only when the divergence of the veloc- 
ity field of the gas surrounding the AGN particle is negative 
- i.e., {V ■ v) < 0. The divergence of the velocity field of the 
neighbour gas particles is calculated using the SPH scheme. 
A constant thermal energy of 10**erg s~^ is deposited to 
the neighbour gas particles when the above condition is sat- 
isfied. This en ergy is roughly consistent with observa tional 
estimates (e.g. iBlanton. Sarazin fc McNamara 1 l2003l) . The 
energy is deposited to the neighbour gas particles in the 
same way as SNe feedback, i.e. smoothed over the neigh- 
bour gas par ticles using the SPH smoo thing algorithm (see 
Sec. 2.3.5 of iKawata fc Gibson ll2003al . for details). Hence, 



^ Thermal energy of 10^'^ erg is applied per supernova, admit- 
tedly unrealistic in terms of energetics, but we are only using the 
model for comparison purposes. 

^ We acknowledge that "turning on" AGN heating at redshift 
z=l "misses" the peak of AGN activity in the Universe by 
~2 Gyr; our choice of z=l for this particular model was not 
driven by empirical cosmological AGN activity arguments, but 
by identifying the redshift for this model at which a stable "nu- 
cleus" could be identified. We will explore this limitation in a 
future study. 
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we assume that the energy of the AGN heating is distributed 
spherically in the scale of the smoothing length (~ kpc) as 
purely thermal energy. This is obviously too simple, and re- 
quired to be improved. However, we adopt this simple model 
for our first study. 

We examine both the resulting X-ray and optical prop- 
erties of the simulation end-products, comparing them quan- 
titatively with observation. The gas particles in our simu- 
lations carry with them knowledge of the density, tempera- 
ture, and abundances of various heavy elements. Using the 
XSPEC vmekal plasma model, we derive the X-ray spec- 
trum for each gas particle, and synthesise them within the 
assumed aperture. We next generate "fake" spectra with the 
response function of the XMM EPN detector, assuming an 
exposure time (40 ks) and target galaxy distance (17 Mpc). 
Finally, our XPSEC fitting provides the X-ray weighted 
temperatures and abundances of various elements. In the 
next section, we compar e the simulation results with the X- 
ray ob servational data in lMatsushita. Ohashi fc Makishima I 
J200d MOMOO). MOMOO used the aperture radius of four 
times the _B-band effective radius. The B-band surface 
brightness profile for our successful model, i.e. Mo del 2, is 
barel y able to be fitted by the Sersic {r^''") law iSersic I 
Il968fl '^. which provides Re,B ~ 11 kpc. Hence we applied a 
projected aperture radius of i? = 45 ^ 4_Rc,s kpc for both 
models. Conversely, the simulated star particles each carry 
their own age and metallicity "tag" , which enables us to gen- 
erate an optical-to-near infrared spectral energy distribution 
for the target galajcy, when combined with our population 
synthesis code (wh ich itself is based upon th e population 
synthesis models of lKodama fc Arimoto l[l997l) . 



3 DISCUSSION AND CONCLUSIONS 

To study both the X-ray and optical properties, we exam- 
ine the Lx — Tx and Lx — [Fe/H]x relations for the X-ray 
properties, and the colour-magnitude relation (CMR) for the 
optical properties, following KG03b. Throughout this paper, 
we norm alise all abundances to the solar "meteoritic" val- 
ues from lAnders fc Grevesse I ||l989) . Figs. and |5| show 
the Lx — Tx and Lx — [Fe/H]x relations and the CMRs for 
Models 1 and 2 at z = 0, compared that observed. Re- 
call, Model 1 corresponds to Model C of KG03b; such mod- 
els are broadly consistent with the empirical X-ray prop- 
erties of ellipticals and agreement driven primarily by ra- 
diative cooling (see alsolFuiita. Fukumoto fc Qkoshi Ill996l: 



diative cooling (see alsoll<uiita. -bukumoto fe Ukoshi 111991 
Pearce et al. II2OO0I: iMuanwong et al. I|200ll: iBorgani et al. 



20021: iD ave. Katz fc Weinberg I 120021: IValdarnini I l2003f 

Tornator e ct al. 2004). 

KG03b showed that radiative cooling ensures that the 
hot dense gas turns into the cold (i.e. non- X-ray-emitting) 
gas, and keeps X-ray-emitting gas high temperature and low 
density. KGOSb also found that stars preferentially enrich 
the gas in the central region where radiative cooling is effi- 
cient. The high-density cold gas is incorporated into future 
generations of stars, and thus a large fraction of the iron 

^ As mentioned in KGOSb, the B-band surface brightness profile 
for Model 1, i.e. Model C in KGOSb, is too much centrally con- 
centrated due to the continuous star formation at the centre, and 
cannot be fitted by the Sersic law. 



> 



□ 



-24 



-22 



>s< X 



-20 



Figure 2. Comparison of the simulated U — V CMR 
(square/circle/triagle for Models 1/2/2L, respectively) and that 
of the C oma cluster ellipticals (cross es). The observational data 
are from lBower. Lucev fc E'iiis^ill992^ . The dashed line shows the 
CMR fitted to the Coma Cluster galaxies. For comparison, the 
solid star co rresponds to NGC 3923 whose total magnitude is ob- 
tained from'de Vaucoulcurs et ahl ll99ll) . and colours are from 
[Pcrsson. Frogel & Aaronsor^~ ^^79l) who presented the colours 
within a aperture similar in size to what iBower. Lucev fc Ellis I 
a99a) used. 



ejected from stars is locked into future generations of stars. 
As a result, the hot gaseous halo which emits X-rays is not 
enriched efficiently, leading to the lower [Fe/H]x for Model 1. 
Hence, radiative cooling ensures that the hot gas has a lower 
metallicity than the stellar component. However, the optical 
colours of the resulting stellar component for this model are 
inconsistent with the observational data (Fig. |5J , due to the 
excessive star formation at low redshift (Fig.|2|l being driven 
by radiative cooling (see KG03b for more details). 

In contrast, the new AGN heating model (Model 2) 
roughly reproduces both the X-ray and optical observational 
data. The shaded area in Fig. displays the range of Lx 
and Tx for a sample of elliptical galaxies with low-Lx and 
high-Tx (group 5 in MOMOO), and 90 % confidence hmit of 
[Fe/H]x from the composite X-ray spectrum of these galax- 
ies. The X-ray properties for Model 2 are consistent with 
those for the sample ellipticals. In addition. Fig. |21 demon- 
strates that the colour and magnitude for Model 2 are con- 
sistent with the CMR of the Coma cluster ellipticals and 
roughly reproduce those for NGC 3923* which is the bright- 
est galaxy in this sample. 

Of greatest importance for the analysis here is the fact 
that this AGN heating is sufficiently efficient to suppress 
the late-time star formation which plagued our earlier study 
(KGOSb). Fig. El demonstrates this graphically. As a result, 
Model 2 succeeds in reproducing - simultaneously - both the 
X-ray and optical properties^ observed in nearby ellipticals 

* NGC 3923 is known as a shell elliptical. However, their X-ray 
propertie s seem to be normal, and old stellar population is domi- 
nant (e.g. lMichard fc Prugniel Il2004) . Note that the effective ra- 
dius of NGC 3923 (ile.s ~6.6 kpc) is smaller than that for Model 
2, which might affect the comparison of the X-ray properties, be- 
cause we adopted the aperture radius of AR^ g . 
^ Fig. 13 shows that Model 2 also reproduces the Type la SNe 
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Figure 3. Time variation of the star formation rate (upper) 
and the event rate of SNe II (thin lines) and SNe la (thick line) 
for Models 1 (left panels), 2 (solid lines in right panels), and 2L 
(dashed lines in right panels). Solid square (circle) with error- bar 
in lowe r left (right) panel is taken from the observational SNe la 
rate bv lCaperrallo. Evans &: Turattol 119991) . taking into account 
the B-band luminosity for Model 1 (2). To show them clearly, 
those SNe la rates are plotted at t = 13 Gyr. 
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Figure 4. Histories of the feedback energy from Type II SNe 
(dashed). Type la SNe (dotted) and the AGN (solid), for Model 2. 

such as NGC 3923. It is also worth noting that the lower 
[Fe/H]x in Model 2 is induced by a different mechanism, 
i.e. blowing out the enriched gas, from that in Model 1, i.e. 
radiative cooling, because the iron ejected from stars is no 
long able to be hidden in the cold gas or future generations 
of stars in Model 2. In fact, we found that compared with 
Model 1, a larger fraction of iron ejected from stars is blown 
out from the system by the AGN heating. 

Fig. 2] shows the histories of the feedback energy from 
SNe II, SNe la and the AGN. Until star formation ceases at 
0~O.4, the heating due to SNe II is comparable to that of 
the AGN. However, once star formation stops, the number 
of SNe II drops, and consequently they become unimportant 

rate observed in nearby ellipticals iCaperrallo. Evans &: Turatto I 
Il999f) within the associated observational uncertainties. 



as a heat source. Conversely, the AGN provides continuous 
heat to the g the AGN activity is not directly linked to 
the star formation well outside the nucleus. SNe la heating 
by itself is too inefficient to prevent the gas from cooling 
and subsequent star formation occurring in Model 1. The 
mean heating energies from SNe la and the AGN over the 
final 3 Gyr of the simulation are 4.8 x 10*^ and 2.5 x lO"^ 
ergs, respectively. The efficient, star formation-independent, 
heating provided by the AGB is required to suppress late- 
time star formation and, consequently, recover the red stellar 
populations of ellipticals. 

Even after star formation ceases at z~0.4, the AGN ac- 
tivity becomes more sporadic, because the gas fueling is also 
suppressed by the AGN heating itself.® This self-regulation 
of the AGN activity ensures that the hot gas present is not 
entirely removed from the system due to over-heating, or 
entirely converted to cold gas by efficient radiative cooling. 
For comparison, we also ran a model in which the AGN 
was continuously active (from redshift z=l to the present- 
day) independent of gas kinematics in the central region. We 
found that such continuous heating blows out the gas from 
the system completely because of the efficient (and uninter- 
rupted) heating. 

Fig. |5] displays the time variation of the relevant X-ray 
properties for our simulated elliptical galaxies. Here, Tx and 
[Fe/H]x are calculated by taking the mean values weighted 
by p^T^/^ for the gas with T > 10*^ K, following Fig. 12 of 
KG03b. Prior to the cessation of star formation (at z~QA), 
the X-ray properties for Model 2 vary sporadically, but settle 
into relevant stability thereafter. We suspect that this may 
be traced to our (too simple) assumption that the AGN be- 
comes active instantaneously at z=l\ recall that this heating 
source was introduced when the simulation was dynamically 
stable and in possession of a deep potential well. As such, 
there was a finite time (^^2-3 Gyr) before the AGN could 
contribute sufficiently to the heat of the surrounding am- 
bient ISM and a new stable configuration reached. In real- 
ity, of course, it is likely that the AGN develops before the 
host system becomes stable, and the AGN heating could 
affect the hot gas more efficiently. In that picture, we can 
anticipate the hot gas properties evolving more passively 
prior at 2>0.4, in contrast with the current models. Unfor- 
tunately, at this time, it is difficult to construct a more self- 
consistent model using three-dimensional cosmological hy- 
drodynamical simulations. Nevertheless, our simple model 
demonstrates that the AGN heating induced by gas infall 
leads to a stable condition for a lengthy period of time due 
to its self-regulated nature. Hence, our results shown in this 
paper are not ephemeral, but represent stable, long-term 
properties of the simulated ellipticals. 

After the hot gas properties reached a stable condition, 
the X-ray properties are still varying a little. This temporal 
variance is considered to represent the range of uncertainties 



® Note that our assumed AGN heating energy is 10** erg s^*. 
Hence, in Fig. |1] AGN heating energy of less than 10** erg s~* 
means that the AGN activity has been "switched on and off" 
within the period of the size of that bin. This is admittedly un- 
realistic, as the bin size is small (typically ~ 40 Myr), compared 
to the typical lifetime of the AGN activity (~100 Myr). This is a 
current limitation of the phenomenological AGN model employed. 
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Figure 5. Time variation of tiie X-ray luminosity (upper left), 
temperature (lower left), iron abundance (upper right) and the 
hot (thin lines) and total baryon (thick line) mass (lower right) 
within R = 45 kpc. The values obtained in Models 1, 2, and 2L 
are presented as dotted, solid, and dashed lines, respectively. 

in our simulation^. Thus, Fig. |3 means that our simulation 
cannot tell any difference in Tx or [Fe/H]x between Models 
1 and 2. Conversely, the difference in Lx is significant. The 
range of Log(Lx) since z = 0.2 {t ki 11 Gyr) is 41.2—41.8 and 
40.0—40.7 for Models 1 and 2, respectively. Consequently, 
the AGN heating model reduces the X-ray luminosity for the 
target galaxy by at least a factor of five, and leads to the X- 
ray properties similar to low-Lx and high-Tx ellipticals . It 
would be interesting to examine whether AGN heating can 
produce the other type of elliptical galaxies, such as high- 
Lx ellipticals, because it seems that there is no difference 
in optical properties between lo w-Lx and high-Lx ellipticals 
jKodama fc MatsushitTlbOOOl) . 

Finally, to demonstrate the uncertainty by the numer- 
ical resolution, we carried out a simulation (Model 2L) in 
which the high-resolution region (Sec. |5J consists of lower- 
resolution particles (27/64 by mass), and the same parame- 
ters as Model 2 are adopted. Note that Model 2L has a dif- 
ferent pattern of density fluctuations in the high-resolution 
region, which leads to a different minor merger history from 
Model 2. Fig. shows that Model 2L provides similar X- 
ray properties to Model 2 within the range of the temporal 
variance seen in Fig. |S1 i.e. this variance corresponds to the 
uncertainties of our simulation as described above. Star for- 
mation is quenched by the AGN heating at 2; ~ 0.3 (Fig. 
ISJ, and the red colour is reproduced (Fig. |2J. Hence, our 
conclusion shown above is less sensitive to the numerical 
resolution. 

Our cosmological simulation with the AGN heating has 
successfully produced a system which is consistent with 
observed elliptical galaxies in both the X-ray and optical 
regimes. To our knowledge, this is the first elliptical galaxy 
model which does so self- consistently. We admit that our 
AGN model is rather simple and our study limited by its 



Some features like a periodic change of Lx in Model 2 seem 
to be real, because it is caused by the expansion and re-collapse 
of the hot gas due to the self-regulated AGN activity. However, 
we postpone such discussion until higher resolution simulations 
become available. 



focus upon only one (representative) elliptical galaxy simu- 
lation. Nevertheless, we believe that this encouraging result 
provides additional impetus for galactic astronomy commu- 
nity to further explore the role that AGN heating plays in 
the formation and evolution of elliptical galaxies. 
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